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(57) Abstract: The invention relates to an imagining method which is used to observe the simultaneous propagation of a shear wave 
^ at a multitude of points in a viscoelastic diffusing medium (1). The inventive method consists in: generating the shear wave by using 

a network of transconductors (6) to emit at least one focused ultrasonic compression wave in the viscoelastic medium; subsequently, 
f^| using the same network of transconductors to emit non-focused ultrasonic compression waves at a rapid rate, in order to produce a 

O scries of images of the medium; and processing the images thus obtained in delayed time in order to determine movements in relation 
tc the :r.cdium during the propagation of the shear wave. 
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(57) Abrege : Proced^ dMmagerie pour observer la propagation d*une onde de cisaillement simultan^ment en une multitude de 
points d'un milieu viscoelastique dififiisant (1). On fait g^n^rer Tonde de cisaillement en tirant au moins une onde ultrasonore de 
compression focalis^e dans le milieu viscoelastique avec un r^seau de transconducteurs (6) , puis on ^met k cadence rapide, avec le 
meme rdseau de transconducteurs, des ondes ultrasonores de compression non focalis^es qui permettent d'obtenir une succession 
d' images du milieu, et on traite en temps diffi^ les images ainsi obtenues pour determiner les mouvements du milieu lors de la 
propagation de Tonde de cisaillement. 
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Procede et dispositif d^imagerie utiliscint des ondes de ci- 
saillemeni:. 



La pr6sente invention est relative aux proc6d6s et 
5 dispositifs d'imagerie utilisant des ondes de cisaillement . 

Plus particuliferement, 1 ' invention concerne un pro- 
cede d'imagerie utilisant des ondes de cisaillement pour 
observer un milieu viscoelastique diffusant qui contient 
des particules ref lechissant les ondes ultrasonores de com- 
10 pression, procede comprenant : 

(a) une etape d' excitation au cours de laquelle on 
genere une onde elastique de cisaillement dans le milieu 
viscoelastique, 

(b) une etape d' observation au cours de laquelle on 
15 observe la propagation de 1 ' onde de cisaillement simultane- 

ment en une multitude de points d'un champ d' observation 
dans le milieu viscoelastique, cette etape d' observation 
comprenant les sous-6tapes suivantes : 

(bl) faire emettre dans le milieu viscoelasti- 
20 que, par un reseau de transducteurs commandes independam- 
ment les uns des autres, une succession de tirs d' ondes ul- 
trasonores de compression non focalis6es a une cadence d'' au 
moins 500 tirs par seconde, 

(b2) faire detecter et enregistrer en temps reel 
25 des signaux acoustiques regus du milieu visco61astique 2, 
comprenant les §chos g6neres par les ondes ultrasonores de 
compression non focalis^es en interagissant avec les parti- 
cules ref lechissantes dudit milieu viscoelastique, 

(c) et au moins une 6tape de traitement au cours de 
30 laquelle : 

(cl) on traite les signaux acoustiques succes- 
sifs regus du milieu viscoelastique au cours de la sous- 
etape (b2) pour determiner des images de propagation suc- 
cessives de I'onde de cisaillement. 
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(c2) et pour determiner au moins un parametre de 
mouvement du milieu visco61astique en differents points du 
champ d' observation. 

On obtient ainsi un "film" illustrant clairement la 

5 propagation de I'onde de cisaillement dans le milieu vis- 
co61astique, qui peut permettre une analyse qualitative 
et/ou quantitative pour reperer des zones de duret6 diffe- 
rente du reste du milieu viscoelastique ou des zones ayant 
un temps de relaxation different du reste du milieu viscoe- 

10 lastique. 

Le document WO-A-00/55 616 decrit un example d'un 
tel precede, dans lequel on genere 1 ' onde de cisaillement 
a la surface du milieu viscoelastique, Ce procede donne 
toute satisfaction notamment pour imager des zones situees 

15 relativement pres de la surface du milieu viscoelastique. 
Mais ce procede connu ne permet pas d' observer certaines 
zones du milieu viscoelastique/ notamment : 

- des zones suffisamment profondes pour ne pas 
pouvoir §tre atteintes par des ondes de cisaillement g§ne- 

20 r6es en surface (les ondes de cisaillement s'attenuent ra- 
pidement) , 

- et des zones d' ombre masqu§es par des obstacles 
(notamment des parties du squelette d'un patient ou des zo- 
nes liquides telles que des kystes liquides) qui gSnent la 

25 propagation des ondes de cisaillement. 

De plus, si le champ d' observation est partielle- 
ment dans une zone d' ombre, il peut etre necessaire de de- 
placer le dispositif g6nerateur- des ondes de cisaillement 
en cours d' observation, ce qui est fastidieux pour I'utili- 

30 sateur. 

Enfin, le dispositif g^nerateur des ondes de ci- 
saillement est relativement lourd et complique le disposi- 
tif. 

La presente invention a notamment pour but de pal- 
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lier ces inconv6nients , 

A cet effet, selon 1' invention, un procede du genre 
en question est caract6ris6 en ce qu'au cours de I'etape 
d' excitation (a), on fait g6nerer I'onde 61astique de ci- 
5 saillement en faisant eitiettre au moins une onde ultrasonore 
focalisee dans le milieu viscoelastique par ledit reseau de 
transducteurs, la focalisation et la chronologie de ladite 
onde ultrasonore focalisee, ainsi que la chronologie desdi- 
tes ondes ultrasonores non focalisees, etant adaptees pour 
10 qu'au moins certaines desdites ondes ultrasonores non foca- 
lisees parviennent dans le champ d' observation lors de la 
propagation de I'onde de cisaillement dans ce champ 
d' observation - 

Ainsi, le meme reseau de transducteurs permet a la 
15 fois d'engendrer I'onde elastique de cisaillement de fagon 
choisie dans le champ d^ observation, et d' observer ensuite 
cette propagation, grSce au fait que le dispositif 
d' imagerie est adapt6 pour g6nerer soit des ondes ultraso- 
nores focalis6es permettant de gen6rer I'onde elastique de 
20 cisaillement, soit des ondes ultrasonores non focalisees 
permettant d' observer la propagation de I'onde de cisaille- 
ment, et grace au choix judicieux : 

de la chronologie de ces differentes Emissions, 
et du ou des points de focalisation de I'onde 
25 ultrasonore focalisee. 

Le procede d' imagerie selon 1' invention est done 
aise k mettre en oeuvre pour un utilisateur, avec un dispo- 
sitif relativement simple et leger. L' invention est done 
d'un tres faible cout de revient en comparaison de techni- 
30 ques concurrentes telles que I'IRM, et elle permet le cas 
echeant de constituer des systemes d' imagerie ambulatoires 
utilisables tant en imagerie preoperatoire qu' en imagerie 
postoperatoire, voire en imagerie per-operatoire - 

A titre d'exemple, dans des applications medicales. 
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le procede selon 1' invention peut permettre de reperer ef- 
ficacement des zones cancereuses dans les tissus d'un pa- 
tient. La propagation des ondes de cisaillement s'y d§roule 
en effet tres dif f 6reinment des zones voisines. Ce reperage 

5 s'effectue beaucoup plus facilement que par une observation 
classique par simple echographie ultrasonore^ puisque la 
propagation des ondes de cisaillement est fonction du mo- 
dule de cisaillement du milieu^ lui-meme tres variable en- 
tre une zone de tissus sains et une zone de tissus cance- 

10 reux : le module de cisaillement varie typiquement dans un 
rapport de 1 a 30 entre une zone saine et une zone cance- 
reuse, alors que le module de compression, qui regit la 
propagation des ondes acoustiques de compression utilisees 
dans 1 ' echographie ultrasonore, varie seulement de I'ordre 

15 de 5% entre un tissu sain et un tissu cancereux. 

De m§me, on peut ainsi reperer des zones necrosees 
au sein d'un tissu, par exemple des zones tumorales ayant 
subi un traitement par hyperthermie aux ultrasons, notam- 
ment pour ^valuer I'efficacite de ce traitement par hyper- 

20 thermie . 

Une autre application possible de 1' invention 
concerne 1' evaluation quantitative du degre de fibrose du 
foie, qui est un parametre important dans les pathologies 
du foie, notamment 1' hepatite C. 

25 On notera que 1' invention permet de generer I'onde 

de cisaillement, et d' observer sa propagation, y compris a 
travers une zone liquide ou une barriere osseuse masquant 
totalement ou partiellement le champ d' observation (crane, 
cage thoracique, etc.)/ puisqu' il est possible de focaliser 

30 des ondes ultrasonores au travers de telles barrieres (voir 
notamment le document WO-A-02/32316 ou la demande de brevet 
frangais n"" 02 10682 du 28 aout 2002) . 

Dans des modes de realisation preferes du procede 
selon 1* invention, on peut eventuellement avoir recours en 
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outre ^ I'une et/ou ^ 1' autre des dispositions suivantes : 

au cours de la sous-etape (b2) , pour determiner 
ledit parametre de mouvement^ on compare plusieurs images 
successives de propagation (par exemple par correlation, 
5 Doppler, etc.) avec une m§me image de reference du milieu 
visco61astique, cette image de reference 6tait d6terminee 
en tirant au moins une onde ultrasonore de compression non 
focalisee dans ledit milieu viscoelastique puis en detec- 
tant et en enregistrant des 6chos generes par ladite onde 

10 ultrasonore de compression non focalisee lorsqu'elle inte- 
ragit avec les particules ref l§chissantes du milieu viscoe- 
lastique (on ameliore ainsi la precision de la mesure du 
parametre de mouvement (par exemple le deplacement) du mi- 
lieu viscoelastique, notamment pour de faibles amplitudes 

15 de mouvement (typiquement moins de 30 \xm avec la technique 
d' excitation par ondes acoustiques utilisee ici) ; 

I'etape (a) est precedee d'une etape d' observa- 
tion initiale (aO) au cours de laquelle on tire au moins 
une onde ultrasonore de compression non focalisee puis on 

20 detecte et on enregistre des echos generes par ladite onde 
ultrasonore de compression non focalisee en interagissant 
avec les particules r6f lechissantes du milieu viscoelasti- 
que, ces echos correspondant (directement ou indirectement ) 
a une image initiale du milieu viscoelastique, et au cours 

25 de la sous-etape (b2) ladite image initiale constitue la- 
dite image de reference pour traiter au moins certaines des 
images successives de deplacement ; 

au cours de I'^tape d' observation initiale (aO) , 
on tire successivement plusieurs ondes ultrasonores de com- 

30 pression non focalisees puis on detecte et on enregistre 
des 6chos generes par chaque onde ultrasonore de compres- 
sion non focalisee en interagissant avec les particules r6- 
f lechissantes du milieu viscoelastique, ces 6chos corres- 
pondant (directement ou indirectement) a plusieurs images 
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pr61iminaires successives du milieu viscoelastique, et on 
determine ladite image initiale du milieu visco61astique en 
combinant lesdites images preliminaires successives ; 

- ledit parametre de mouvement est un deplacement 
5 du milieu visco61astique ; 

I'onde ultrasonore focalis6e 6mise au cours de 
I'etape d' excitation (a) presente une frequence f comprise 
entre 0,5 et 15 MHz et est emise pendant une duree k/f en 
secondes, ou k est un entier compris entre 50 et 5000 et f 
10 est exprime en Hz ; 

I'onde ultrasonore focalisee emise au cours de 
I'etape d' excitation (a) presente une frequence comprise 
entre 0,5 et 15 MHz et est emise pendant une succession de 
periodes d' emission separees par des periodes de repos, les 
15 periodes d' emission se succedant a une cadence comprise en- 
tre 10 et 1000 emissions par seconde ; 

I'onde ultrasonore focalisee emise au cours de 
I'etape d' excitation (a) est une combinaison lineaire (no- 
tarament une somme) de deux signaux monochromatiques de fre- 
20 quences respectives fl et f2 telles que 20 Hz ^ | f 1 - f 2 | ^ 
1000 Hz ; 

I'onde ultrasonore focalis6e 6mise au cours de 
I'etape d' excitation (a) est focalisee simultanement en 
plusieurs points ; 

25 - I'etape (c) de traitement d' images est suivie 

(imm^diatement ou non) par une etape (d) de cartographie au 
cours de laquelle, a partir d'une evolution du parametre de 
mouvement au cours du temps, on calcule au moins un parame- 
tre de propagation de I'onde de cisaillement en au moins 

30 certains points du champ d' observation pour determiner ain- 
si une cartographie dudit parametre de propagation dans le 
champ d' observation ; 

- le parametre de propagation de I'onde de cisail- 
lement qui est calcule au cours de 1' etape de cartographie 
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est choisi parmi la vitesse des ondes de cisaillement, le 
module de cisaillement, - le module d'Young^ 1' attenuation 
des ondes de cisaillement, 1' elasticity de cisaillement^ la 
viscosite de cisaillement et le temps de relaxation mecani- 
5 que ; 

on repete successivement les etapes (a) a (d) en 
emettant des ondes ultrasonores differentes au cours des 
6tapes d' excitation (a) successives, puis on combine les 
differentes cartographies obtenues au cours des etapes de 

10 cartographie (d) successives, pour calculer une cartogra- 
phie synthetique du champ d' observation . 

Par ailleurs, 1' invention a egalement pour objet un 
dispositif d'imagerie utilisant des ondes de cisaillement 
pour observer un milieu viscoelastique diffusant qui 

15 . contient des particules ref lechissant les ondes ultrasono- 
res de compression, ce dispositif comprenant un reseau de 
transducteurs commandes independamment les uns des autres 
par au moins une unite centrale electronique adaptee pour : 

- f aire generer au moins une onde elastique de ci- 
20 saillement dans le milieu viscoelastique, 

observer la propagation de I'onde de cisaille- 
ment simultan6ment en une multitude de points d'un champ 
d' observation dans le milieu viscoelastique, en faisant 
emettre dans le milieu viscoelastique, par ledit reseau de 

25 transducteurs, une succession de tirs d' ondes ultrasonores 
de compression non focalisees, ci une cadence d' au moins 500 
tirs par seconde, puis en faisant dfetecter en temps reel 
par ledit r§seau de transducteurs, et pour enregistrer en 
temps r6el, des signaux acoustiques regus du milieu viscoe- 

30 lastique, comprenant les 6chos g6ner§s par les ondes ultra- 
sonores de compression non focalisees en interagissant avec 
les particules ref 16chissantes dudit milieu visco61astique, 

- et traiter les signaux acoustiques successifs 
re?us du milieu viscoelastique pour determiner des images 
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de propagation successives de I'onde de cisaillement, puis 
determiner au moins un parametre de mouvement du milieu 
viscoelastique en differents points du champ d' observation, 
caract6rise en ce que 1' unite centrale 61ectronique est 

5 adaptee pour faire gfenerer I'onde §lastique de cisaillement 
en faisant emettre au moins une onde ultrasonore focalis6e 
dans le milieu viscoelastique par ledit r6seau de transduc- 
teurs, la focalisation et la chronologie de ladite onde ul- 
trasonore focalis6e, ainsi que la chronologie desdites on- 

10 des ultrasonores non focalis6es, etant adaptees pour que 
lesdites ondes ultrasonores non focalisees parviennent dans 
le champ d' observation lors de la propagation de I'onde de 
cisaillement dans ce champ d' observation . 

D'autres caracteristiques et avantages de 1 ' inven- 

15 tion apparaltront au cours de la description suivante d'une 
de ses formes de realisation, donnee a titre d'exemple non 
limitatif, en regard du dessin joint, 

Sur le dessin, la figure 1 est une vue schematique 
d'un dispositif d'imagerie par ondes de cisaillement selon 

20 une forme de realisation de 1' invention. 

Le dispositif d'imagerie 1 represents sur la figure 
1 est destine a etudier la propagation des ondes 61astiques 
de cisaillement dans un milieu viscoelastique 2 qui est 
diffusant vis ^ vis des ondes ultrasonores de compression, 

25 et qui peut §tre par exemple : 

un corps inerte, notamment dans le cas du 
controle de quality pour des applications industrielles, 

- ou un corps vivant, par exemple une partie du 
corps d'un patient, dans le cas des .applications m6dicales. 

30 Ces mouvement s sont suivis par exemple au mo yen 

d'un micro-ordinateur 4 (comprenant au moins une interface 
d' entree 4a telle qu'un clavier ou autre, et une interface 
de sortie telle qu'un ecran ou autre) ou toute autre unit 6 
centrale electronique, qui fait envoyer dans le milieu 2, a 
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partir de sa surface exterieure 3, des ondes ultrasonores 
de compression qui interagissent avec les particules diffu- 
santes 5 contenues dans le milieu 1/ lesquelles particules 
sont r§f lechissantes pour les ondes ultrasonores de com- 
5 pression. Les particules 5 peuvent §tre constitutes par 
toute heterogeneite du milieu 1, et notamment, lorsqu^il 
s'agit d'une application medicale, par des particules de 
collagene presentes dans les tissus humains (ces particules 
forment sur les images echographiques des points connus 

10 sous le terme "speckle") . 

Pour observer la propagation de I'onde de cisaille- 
ment^ on utilise une sonde ultrasonore 6 disposee centre la 
surface exterieure 3 du milieu observe 1. Cette sonde en- 
voie, selon un axe des impulsions d' ondes ultrasonores 

15 de compression du type de celles couramment utilisees en 
echographies a une frequence comprise par exemple entre 0,5 
et 100 MHz et de preference entre 0,5 et 15 MHz, par exem- 
ple de 1 ' ordre de 4 MHz. 

La sonde ultrasonore 6 est constitute par un reseau 

20 de n transducteurs ultrasonores Tl, T2, Ti, - - . Tn, n 

etant un nombre entier n au moins egal a 1. 

Cette sonde 6 peut se presenter par exemple sous la 
forme d'une barrette lineaire pouvant comprendre par exem- 
ple n = 128 transducteurs alignes selon un axe Y perpendi- 

25 culaire a I'axe X. Mais la sonde en question pourrait §tre 
egalement un reseau bidimensionnel (plan ou non) de trans- 
ducteurs . 

Les transducteurs Tl, T2, ... Tn sont commandes in- 
dependamment les uns des autres par le micro-ordinateur 4, 
30 §ventuellement par 1 ' interm6diaire d'une unite centrale CPU 
qui est contenue par exemple dans une bale electronique 7 
relive par un cSble souple a la sonde 6. Les transducteurs 
Tl-Tn peuvent ainsi emettre selectivement : 

- soit une onde ultrasonore de compression "plane" 
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(c'est k dire en 1' occurrence une onde dont le front d'onde 
est rectiligne dans le plan X, Y) ou tout autre type d'onde 
non focalisfee eclairant 1' ensemble du champ d' observation 
dans le milieu 2, par exemple une onde generee en faisant 
5 emettre des signaux acoustiques aleatoires par les diffe- 
rents transducteurs Tl-Tn^ 

- soit une onde ultrasonore de compression focali- 
see en un ou plusieurs points du milieu 2. 

Pour observer la propagation de I'onde de cisaille- 
10 ment dans le milieu 2, on procede en plusieurs etapes sue- 
cessives : 

(a) une etape d' excitation au cours de laquelle le 
micro-ordinateur 4 fait generer une onde elastique de ci- 
saillement dans le milieu viscoelastique 2, en faisant 

15 emettre au moins une onde ultrasonore focalisee dans le mi- 
lieu viscoelastique par la sonde 6, 

(b) une etape d' observation au cours de laquelle on 
observe la propagation de 1 ' onde de cisaillement simultane- 
ment en une multitude de points du champ d * observation dans 

20 le milieu viscoelastique 2, cette etape d' observation com- 
prenant les sous-etapes suivantes : 

(bl) le micro-ordinateur 4 fait emettre dans 
le milieu viscoelastique, par la sonde 6 une succession de 
tirs d'ondes ultrasonores de compression non focalis6es a 

25 une cadence d' au moins 500 tirs par seconde (la focalisa- 
tion et la chronologie de l*onde ultrasonore focalis6e 
6mise a 1* etape (a), ainsi que la chronologie desdites on- 
des ultrasonores non focalis6es, sont adaptees pour qu'au 
moins certaines desdites ondes ultrasonores non focalisees 

30 parviennent dans le champ d' observation lors de la propaga- 
tion de I'onde de cisaillement dans ce champ d' observation, 
pour au moins certaines emissions d'onde ultrasonore non 
focalis6e) , 

(b2) le micro-ordinateur 4 fait detecter par la 
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sonde 6 et enregistrer en temps reel des signaux acousti- 
ques re<?us du milieu viscoelastique 2, comprenant les echos 
gen6r6s par les ondes ultrasonores de compression non foca- 
lis6es en interagissant avec les particules r6f 16chissantes 
5 5 du milieu viscoelastique, ces echos correspondant (direc- 
tement ou indirectement) a des images successives de depla- 
cement du milieu viscoelastique, 

(c) et au moins une etape de traitement au cours de 
laquelle : 

10 (cl) le micro-ordinateur 4 traite les signaux 

acoustiques successifs regus du milieu viscoelastique 2 au 
cours de la sous-etape (b2) pour determiner des images de 
propagation successives, 

(c2) et le micro-ordinateur 4 determine au moins 

15 un parametre de mouvement du milieu viscoelastique 2 en 
differents points du champ d' observation . 

L'onde ultrasonore focalisee emise au cours de 
1' etape d' excitation (a) peut etre une onde monochromatique 
de frequence f comprise entre 0,5 et 15 MHz, par exemple 

20 environ 4 MHz, emise pendant une duree k/f en secondes, ou 
k est un entier compris entre 50 et 5000 (par exemple de 
I'ordre de 500) et f est exprim^ en Hz. Eventuellement , une 
telle onde peut §tre emise pendant une succession de perio- 
des d' emission separees par des periodes de repos, les pe- 

25 riodes d' Emission se succedant a une cadence comprise entre 
10 et 1000 Emissions par seconde. 

En variante, l^onde ultrasonore focalisee emise au 
cours de 1' etape d' excitation (a) est une combinaison li- 
neaire (notamment une somme) de deux signaux monochromati- 

30 ques de frequences respectives fl et f2 telles que 20 Hz ^ 
|fl - f2| ^ 1000 Hz, ce qui produit une onde modul6e en am- 
plitude avec une frequence de modulation | f 1 - f 2 | . 

Par ailleurs, I'onde ultrasonore focalisee emise au 
cours de 1' etape d' excitation (a) peut eventuellement etre 
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focalis6e, simultanement ou non, en plusieurs points^ de 
fagon que I'onde de cisaillement generee presente la forme 
souhait^e (par exemple, on peut ainsi gen6rer une onde de 
cisaillement plane ^ ou au contraire une onde de cisaille- 
5 ment focalis6e) et illumine les zones souhait^es du milieu 
2. 

Au cours de I'etape (bl), qui peut durer par exem- 
ple moins d*une seconde, on peut 6mettre les ondes ultraso- 
nores de compression non focalisees^ a une cadence comprise 

10 entre 500 et 10 000 tirs par seconde et de pr§f6rence com- 
prise entre 1000 et 5 000 tirs par seconde (cette cadence 
est limitee par le temps d ' aller-retour de I'onde de com- 
pression dans le milieu 2, done par I'epaisseur du milieu 2 
dans la direction X : il faut en effet que tous les echos 

15 generes par I'onde de compression aient ete regus par la 
sonde 6 avant d'envoyer une nouvelle onde de compression). 

Chaque onde ultrasonore de compression non focali- 
see se propage dans le milieu 2 avec une vitesse de propa- 
gation beaucoup plus elevee que les ondes de cisaillement 

20 (par example de 1 ' ordre de 1500 m/s dans le corps humain) ^ 
et interagit avec les particules ref lechissantes 5, ce qui 
genere des echos ou autres perturbations analogues du si- 
gnal, connus en soi sous le nom de "bruits de speckle" dans 
le domaine de 1 * echographie . 

25 Ces "bruits de speckle" sont capt6s par les trans- 

ducteurs Tl, ... , Tn au cours de la sous-etape (b2), apres 
chaque tir d'onde ultrasonore de compression non focalisee. 
Le signal sij (t) ainsi capt6 par chaque transducteur Ti 
apres le tir n** j est tout d'abord echantillonne ^ haute 

30 frequence (par exemple de 30 a 100 MHz) et numerise en 
temps r6el (par exemple sur 12 bits) par un 6chantillonneur 
appartenant ^ la bale 7 et reli6 a ce transducteur, respec- 
tivement El, E2, ... En. 

Le signal sij (t) ainsi echantillonne et numerise 
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est ensuite memorise, §galement en temps r6el, dans une m6- 
moire Mi appartenant ^ la baie 7 et propre au transducteur 
Ti. 

Chaque m§moire Mi presente par exemple une capacite 

5 de I'ordre de 128 Mo, et contient 1' ensemble des signaux 
sij (t) re(?us successivement pour les tirs j = 1 a p. 

En temps differe, apres la memorisation de tous les 
signaux sij (t) correspondant a une meme propagation d'onde 
de cisaillement, 1' unite centrale CPU fait retraiter ces 

10 signaux par un circuit sommateur S appartenant a la baie 7 
(ou bien elle effectue elle-meme ce traitement, ou encore 
ledit traitement peut etre effectue dans le micro- 
ordinateur 4), par un processus classique de formation de 
voies correspondant a la sous-etape (cl) . 

15 On genere ainsi des signaux Sj (x, y) qui corres- 

pondent chacun a 1' image du champ d' observation apres le 
tir n** j (dans le cas ou I'onde ultrasonore non focalisee 
est une onde plane) . 

Par exemple, on peut determiner un signal Sj (t) 

20 par la formule suivante : 

SjV) = ^ a, y).sij[t{x, y) + d, (x, y) I V] 

OU : 

- sij est le signal brut pergu par le transducteur 
n® i apres le tir d'onde ultrasonore de compression n° j, 
25 - t(x,y) est le temps mis par I'onde ultrasonore 

de compression pour atteindre le point du champ 
d' observation de coordonn6es (x,y), avec t = 0 au d6but du 
tir n** j, 

di(x,y) est la distance entre le point du champ 
30 d' observation de coordonnees (x,y) et le transducteur n"* i, 
ou une approximation de cette distance, 

V est la Vitesse moyenne de propagation des on- 
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des acoustiques ultrasonores de compression dans le milieu 

viscoelastique observe. 

et ai (X, y) est un coefficient de ponderation 

tenant compte de lois d' apodisation (en pratique, on pourra 
5 dans de nombreux cas considerer que ai(x,y)=l)- 

La formule ci-dessus s' applique mutatis mutandis 

lorsque le champ d' observation est a 3 dimensions (reseau 

bidimensionnel de transducteurs) , en remplagant les coor- 

donnees spatiales (x,y) par (x,y,z). 
10 Apres 1 ' etape eventuelle de formation de voies, 

1' unite centrale CPU memorise dans una memoire centrale M 

appartenant a la bale 7 les signaux d' images Sj{x,y) ou 

Sj (x) ou Sj(x,y,z), qui correspondent chacun au tir n° j. 

Ces signaux peuvent egalement etre memorises dans le micro- 
15 ordinateur 4 lorsqu' il ef f ectue lui-meme le traitement 

d' image. 

Ces images sont ensuite traitees en temps differe a 
la sous-6tape (c2), par correlation et avantageusement par 
intercorr61ation soit deux a deux, soit de preference avec 

20 une image de r6f6rence qui peut etre : 

- soit une image de d^placement precedemment de- 
termin6e comme expliqu6 ci-dessus et utilisee corame image 
de reference pour les images de d6placement ulterieures (ou 
pour un nombre limite d' images de deplacement ulterieures, 

25 par exemple 30 images de deplacement), 

soit d6terminee au cours d'une etape pr§limi- 
naire d' observation initiale (aO) , comme les images succes- 
sives de deplacement susmentionnees, en faisant emettre une 
ou plusieurs ondes ultrasonores non focalisees par la sonde 

30 6 avant 1* etape d' excitation (a) qui g6n6re I'onde de ci- 
saillement (lorsque plusieurs ondes ultrasonores de com- 
pression non focalisees sont ainsi emises avant la phase 
d' excitation, on enregistre des echos generes par chaque 
onde ultrasonore de compression non focalisee en interagis- 
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sant avec les particules ref 16chissantes du milieu viscoe- 
lastique, ces echos correspondant ^ plusieurs images preli- 
minaires successives du milieu viscoelastique, at on deter- 
mine ladite image initiale du milieu viscoelastique en com- 
5 binant lesdites images preliminaires successives et notam- 
ment en moyennant les valeurs des pixels desdites images 
preliminaires) . 

L'intercorrelation susmentionn6e peut §tre r^alisee 
par exemple dans un circuit electronique specialise DSP ap- 

10 partenant a la bale 7, ou etre programmee dans 1' unite cen-- 
trale CPU ou dans le micro-ordinateur 4. 

Au cours de ce processus d' intercorrelation, on 
maximise una fonction d ' intercorrelation 

<Sj (x, y) , Sj+1 (x, y) > afin de determiner le deplacement subi 

15 par chaque particule 5 donnant lieu a un echo ultrasonore. 

Des exemples de tels calculs d' intercorrelation 
sont donnes dans I'etat de la technique^ notainment par 
O'Donnell et al • ("Internal displacement and strain imaging 
using speckle tracking"., IEEE transactions on ultrasonic, 

20 f erroelectrics, and frequency control, vol. 41, n*" 3, mai 
1994, p. 314-325) et par Ophir et al, ("Elastography : a 
quantitative method for imaging the elasticity of biologi- 
cal tissues". Ultrasonic imag., vol, 13, p. 111-134, 1991). 

On obtient ainsi un ensemble de vecteurs deplace- 

25 ments "('^'O engendres par les ondes de cisaillement en cha- 
que position r du milieu 2 sous I'effet de I'onde de ci- 
saillement (ces vecteurs deplacement s peuvent 6ventuelle- 
ment §tre r6duits ^ une seule composante dans 1' exemple 
considere ici) . 

30 Get ensemble de vecteurs d§placements est stocke 

dans la memoire M ou dans le micro-ordinateur 4 et peut par 
exemple etre visualise, notamment au moyen de l'§cran 4a du 
micro-ordinateur, sous la forme d'un film ralenti ou la va- 
leur des deplacements est illustree par un parametre opti- 
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que tel que par un niveau de gris ou par un niveau chroma- 
tique . 

On visualise ainsi parf aitement les differences de 
propagation de I'onde de cisaillement entre les zones de 
5 caracteristiques diff6rentes du milieu 2, par exemple les 
tissus sains et les tissus canc6reux dans le cas d'une 
application inedicale. 

Ce film de propagation de I'onde de cisaillement 
est en outre superposable avec une image echographique 

10 classique, qui peut etre generee egalement par le disposi- 
tif 1 decrit ci-dessus. 

Par ailleurs;- il est egalement possible de calculer 
non pas les deplacements de chaque point du milieu observe 
2, mais les deformations du milieu 2, c'est a dire des vec- 

15 teurs dont les composantes sont les derivees des composan- 
tes des vecteurs deplacements respectivement par rapport 
aux variables d'espace (coordonnees selon X, Y dans 1' exem- 
ple considere) . Ces vecteurs de deformation sont utilisa- 
bles comme les vecteurs deplacements pour visualiser clai- 

20 rement la propagation de I'onde de cisaillement sous la 
forme d'un film, et pr^sentent en outre I'avantage de s'af- 
franchir des deplacements de la sonde 6 par rapport au mi- 
lieu observe 2. 

A partir des champs de deplacements ou de deforma- 

25 tions, le micro-ordinateur 4 peut avantageusement proceder 
ensuite a une etape de cartographie (d) au cours de la- 
quelle, k partir de 1' evolution du parametre de mouvement 
(deplacement ou deformation) au cours du temps dans le 
champ d' observation X, Y (ou X, Y, Z dans le cas d'un re- 

30 seau bidimensionnel de transducteurs) , on calcule au moins 
un parametre de propagation de I'onde de cisaillement, soit 
en certains points du champ d' observation choisis par 
1 'utilisateur a partir du micro-ordinateur 4, soit dans 
tout le champ d' observation. 



wo 2004/021038 



17 



'CT/FR2003/002S16 



Le paraitietre de propagation de I'onde de cisaille- 
ment qui est calcule au cours- .de I'etape de cartographie 
est choisi par exemple parmi : la .vitesse Cs des ondes de 
cisaillement, le module de cisaillement ou le module 
5 d' Young E=3^/ 1' attenuation a des ondes de cisaillement, 
l'61asticite de cisaillement ..jil, la viscosite de cisaille- 
ment ii2, ou le temps de relaxation mecanique xs des- tissiis. 

Par exemple, on peut calculer en diff6rents points 
du champ d' observation : 
10 - la valeur de la c616rite Cs de I'onde de cisail- 

lement,- qui donne accfes k la duret6 des tissus, 

- la* valeur du temps- de relaxation mecanique xs 
des tissus, ' caracteri-stique de la viscosity locale du mi-, 
lieu- 

15 Pour. cela, on utilise I'^quation de propagation (1) 

suivante, ^ laquelle obeissent les deplacements "('"'^) en- 
gendres par les ondes de cisaillement en chaque position ^ 
du milieu : 

20 . - (1) 

ou p est la densite des tissus, xS est le temps de relaxa- 
tion mecanique des tissus et cS est la celerite de I'onde 
de cisaillement, directement. relive au module d'Young E des 
tissus par la relation : 

25 



30 




(2) 



Dans le domaine de Fourier, 1' Equation d'onde (1) 
ci-dessus peut s'ecrire : 

a}^pUir,a>) = c/ (1 AC/(r,ffl) ^ 3 ^ 



10 
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est la transform^ de Fourier du champ de d6place- 
ment "('^'^) mesure en chaque point et ^^('^'^^est la trans- 
f ormee de Fourier du laplacien spatial de ce champ ^^^^^ . 
Etant donne que coxS « 1, on a done acces ^ une expression 
simplifiee : 

.2 U(r,co) 
Co =£y p 

' '^A£/(r,<B) (4) 
<D 6XJ{r,io) (5) 



oii "^C^) est la phase de la variable complexe x. 
Les fonctions ^('^'^^et ^U(f,a}) gtant connues en chaque 
point de 1' image echographique, il est done possible de 
mesurer en tout point de I'espace le module d' Young et le 

15 temps de relaxation mecanique des tissus, en etablissant 
ainsi une cartographie de ces deux parametres. 

De plus^ les equations (4) et (5) etant vraies a 
chaque frequence, le calcul de cS et xS peut etre avanta- 
geusement moyenne sur 1' ensemble de la bande des frequences 

20 portees par I'onde de cisaillement, ameliorant ainsi forte- 
ment la qualite de la cartographie realisee. A cet effet, 
on peut utiliser les formules suivantes : 



Co = co^p — ^ do 



ou coO et CO 1 sont les frequences minimales et maximales que 
30 porte I'onde de cisaillement 
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Le mode de calcul serait le meme en utilisant non 
plus les d6placements, mais les def oarmations du milieu ob- 
serve 2- 

Par ailleurSf on peut avantageusement 6tablir sue- 
5 cessivement plusieurs cartographies des parametres de pro- 
pagation souhait6s, par exemple cS et xS, en generant suc- 
cessivement des ondes de cisaillement differentes, obtenues 
par exemple en 6mettant des ondes ultrasonores de compres- 
sion focalisees successivement en plusieurs points ou ayant 
10 des formes d' ondes differentes. On peut ensuite combiner 
les differentes cartographies obtenues, par exemple en les 
moyennant, de fagon a obtenir une cartographie synthetique 
plus* riche et plus precise. 
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REVEiroiCATIONS 

1. Precede d'imagerie utilisant des ondes de ci- 
saillement pour observer un milieu visco61astique diffusant 
5 (2) qui contient des particules (5) ref lechissant les ondes 
ultrasonores de compression, procede comprenant : 

(a) una etape d' excitation au cours de laquelle on 
genere une onde elastique de cisaillement dans le milieu 
viscoelastique (2), 

10 (b) une etape d' observation au cours de laquelle on 

observe la propagation de I'onde de cisaillement simultane- 
ment en une multitude de points d'un champ d' observation 
dans le milieu viscoelastique (2), cette etape d'observa- 
tion comprenant les sous-6tapes suivantes : 

15 (bl) faire 6mettre dans le milieu viscoelas- 

tique {2) f par un reseau de transducteurs (6) commandes in- 
d6pendamment les uns des autres, une succession de tirs 
d' ondes ultrasonores de compression non focalisees a une 
cadence d' au moins 500 tirs par seconde, 

20 (b2) faire detect er et enregistrer en temps reel 

des signaux acoustiques regus du milieu viscoelastique (2) , 
comprenant les 6chos gen6r6s par les ondes ultrasonores de 
compression non focalisees en interagissant avec les parti- 
cules ref lechissantes (5) dudit milieu viscoelastique/ 

25 (c) et au moins une etape de traitement au cours de 

laquelle : 

(cl) on traite les signaux acoustiques succes- 
sifs regus du milieu viscoelastique (2) au cours de la 
sous-etape {b2) pour determiner des images de propagation 
30 successives de 1 ' onde de cisaillement, 

(c2) et on determine au moins un parametre de 
mouvement du milieu viscoelastique (2) en differents points 
du champ d^ observation, 

caracterise en ce qu' au cours de 1' etape d' excitation (a). 
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on fait generer I'onde elastique de cisaillement en faisant 
emettre au moins une onde ultrasonore focalis6e dans le mi- 
lieu viscoelastique (2) par ledit reseau de transducteurs 
(6) , la focalisation et la chronologie de ladite onde ul- 

5 trasonore focalisee, ainsi que la chronologie desdites on- 
des ultrasonores non focalisees, etant adaptees pour qu'au 
moins certaines desdites ondes ultrasonores non focalis6es 
parviennent dans le champ d' observation lors de la propaga- 
tion de I'onde de cisaillement dans ce champ d' observation, 

10 pour au moins certaines Emissions d'onde ultrasonore non 
f ocalisee . 

2. Precede selon la revendication 1, dans lequel, 
au cours de la sous-etape (b2) , pour determiner ledit para- 
metre de mouvement, on compare plusieurs images successives 

15 de propagation avec une meme image de reference du milieu 
viscoelastique (2) , cette image de reference etait determi- 
nee en tirant au moins une onde ultrasonore de compression 
non focalisee dans ledit milieu viscoelastique puis en de- 
tectant et en enregistrant des echos generes par ladite 

20 onde ultrasonore de compression non focalisee lorsqu'elle 
interagit avec les particules ref lechissantes (5) du milieu 
viscoelastique • 

3. Precede selon la revendication 2, dans lequel 
I'etape (a) est prec6dee d'une etape d' observation initiale 

25 (aO) au cours de laquelle on tire au moins une onde ultra- 
sonore de compression non focalisee puis on detecte et on 
enregistre des 6chos g6n6res par ladite onde ultrasonore de 
compression non focalisee en interagissant avec les parti- 
cules ref lechissantes (5) du milieu viscoelastique^ ces 

30 6chos correspondant k une image initiale du milieu viscoe- 
lastique, et au cours de la sous-etape (b2) ladite image 
initiale constitue ladite image de reference . pour traiter 
au moins certaines des images successives de deplacement. 

4. Precede selon la revendication 3, dans lequel. 
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au cours de l'6tape d' observation initials (aO), on tire 
successivement plusieurs ondes ultrasonores de compression 
non focalisees puis on d6tecte et on enregistre des 6chos 
g6n§res par chaque onde ultrasonore de compression non fo- 
5 calis^e en interagissant avec les particules ref 16chissan- 
tes (5) du milieu viscoelastique, ces echos correspondant k 
plusieurs images successives du milieu viscoelastique, et 
on determine ladite image initiale du milieu visco61astique 
en combinant lesdites images successives. 

10 5. Precede selon I'une quelconque des revendica- 

tions precedentes, dans lequel ledit parametre de mouvement 
est un deplacement du milieu viscoelastique (2) . 

6- Precede selon I'une quelconque des revendica- 
tions precedentes, dans lequel l^onde ultrasonore focalisee 

15 emise au cours de I'etape d' excitation (a) presente une 
frequence f comprise entre 0,5 et 15 MHz et est emise pen- 
dant une duree k/f en secondes, ou k est un entier compris 
entre 50 et 5000 et f est exprime en Hz. 

7. Proc6de selon I'une quelconque des revendica- 
20 tions 1 a 5, dans lequel I'onde ultrasonore focalisee §mise 

au cours de l'6tape d' excitation (a) presente une frequence 
comprise entre 0,5 et 15 MHz et est emise pendant une suc- 
cession de p^riodes d' Emission s6parees par des periodes de 
repos, les periodes d' Emission se succedant ^ une cadence 
25 comprise entre 10 et 1000 emissions par seconde. 

8 . Precede selon 1 ' une quelconque des revendica- 
tions 1 a 5, dans lequel I'onde ultrasonore focalis6e emise 
au cours de I'etape d' excitation (a) est une corabinaison 
lin6aire (notamment une somme) de deux signaux monochroma- 

30 tiques de frequences respectives fl et f2 telles que 20 Hz 
^ |fl - f2| ^ 1000 Hz. 

9. Precede selon I'une quelconque des revendica- 
tions precedentes, dans lequel I'onde ultrasonore focalisee 
emise au cours de I'etape d' excitation (a) est focalisee 
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simultan6ment en plusieurs points. 

10. Precede selon I'une quelconque des revendica- 
tions precedentes/ dans lequel I'etape (c) de traitement 
d' images est suivie par une 6tape (d) de cartographie au 

5 cours de laquelle, k partir d'une evolution du parametre de 
mouvement au cours du temps, on calcule au moins un parame- 
tre de propagation de I'onde de cisaillement en au moins 
certains points du champ d' observation pour determiner ain- 
si une cartographie dudit parametre de propagation dans le 
10 champ d' observation. 

11. Precede selon I'une quelconque des revendica- 
tions pr§c6dentes, dans lequel le parametre de propagation 
de I'onde de cisaillement qui est calcule "au cours de 
I'etape de cartographie (d) est choisi parmi la vitesse des 

15 ondes de cisaillement, le module de cisaillement, le module 
d' Young, 1' attenuation des ondes de cisaillement, 
1' elasticity de cisaillement, la viscosite de cisaillement 
et le temps de relaxation mecanique, 

12. Precede selon la revendication 11, dans lequel 
20 on repete successivement les etapes (a) a (d) en emettant 

des ondes ultrasonores focalisees differentes au cours des 
etapes d' excitation (a) successives, puis on combine les 
differentes cartographies obtenues au cours des etapes de 
cartographie (d) successives, pour calculer une cartogra- 
25 phie synth6tique du champ d' observation . 

13. Dispositif d'imagerie pour la mise en ceuvre 
d'un precede selon I'une quelconque des revendications pre- 
cedentes, utilisant des ondes de cisaillement pour observer 
un milieu visco61astique diffusant (2) qui contient des 

30 particules (5) r6f l^chissant les ondes ultrasonores de com- 
pression, ce dispositif comprenant un r^seau de transduc- 
teurs (6) commandos independamment les uns des autres par 
au moins une unite centrale 61ectrenique (4, CPU) adapt^e 
pour : 
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faire gen6rer au moins une onde 61astique de ci- 
saillement dans le milieu viscoelastique (2), 

observer la propagation de I'onde de cisaille- 
ment simultanement en une multitude de points d'un champ 

5 d' observation dans le milieu viscoelastique (2), en faisant 
emettre dans le milieu viscoelastique, par ledit reseau de 
transducteurs (6), une succession de tirs d'ondes ultraso- 
nores de compression non focalisees, a une cadence d'au 
moins 500 tirs par seconde, puis en faisant detecter en 

10 temps reel par ledit reseau de transducteurs (6), et pour 
enregistrer en temps reel, des signaux acoustiques regus du 
milieu viscoelastique (2), comprenant les echos generes par 
les ondes ultrasonores de compression non focalisees en in- 
teragissant avec les particules ref lechissantes (5) dudit 

15 milieu viscoelastique, 

et traiter les signaux acoustiques successifs 
regus du milieu viscoelastique (2) pour determiner des ima- 
ges de propagation successives de I'onde de cisaillement , 
puis determiner au moins un param^tre de mouvement du mi- 

20 lieu viscoelastique (2) en diff brents points du champ 
d' observation, 

caraGi:erise en ce cjue 1' unite centrale electronique (4) est 
adaptee pour faire generer I'onde 61astique de cisaillement 
en faisant 6mettre au moins une onde ultrasonore focalisee 

25 dans le milieu viscoelastique par ledit reseau de transduc- 
teurs (6), la focalisation et la chronologie de ladite onde 
ultrasonore focalisee, ainsi que la chronologie desdites 
ondes ultrasonores non focalisees, etant adaptees pour que 
lesdites ondes ultrasonores non focalisees parviennent dans 

30 le champ d' observation lors de la propagation de I'onde de 
cisaillement dans ce champ d' observation • 
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